TIFR/TH/01-49 
NUB-TH/3225 



Super symmetric Dark Matter and Yukawa 

o . 

^ Unification 
o 

Q 



Utpal Chattopadhyay^") , Achille Corsetti^'') and Pran Nath^'') 

^""^ Department of Theoretical Physics, 
Tata Institute of Fundamental Research, Homi Bhabha Road 
Mumbai 400005, India 

2 I ''^''Department of Physics, Northeastern University, Boston, MA 02115, USA 

O 

1^ ' Abstract 
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An analysis of supersymmetric dark matter under the Yukawa unification constraint is 
D , given. The analysis utilizes the recently discovered region of the parameter space of models 



> 

X 



with gaugino mass nonuniversalities where large negative supersymmetric corrections to 
the b quark mass appear to allow h — r unification for a positive sign consistent with the 
^ i 6 s + 7 and (7^ — 2 constraints. In the present analysis we use the revised theoretical 

determination of a^^^^ {a^ = {g^ — 2)/2) in computing the difference a^^^ — a^^^ which takes 
account of a reevaluation of the light by light contribution which has a positive sign. The 
analysis shows that the region of the parameter space with nonuniversalities of the gaugino 
masses which allows for unification of Yukawa couplings also contains regions which allow 
satisfaction of the relic density constraint. Specifically we find that the lightest neutralino 
mass consistent with the relic density constraint, br unification for SU(5) and b — t — t 
unification for SO(IO) in addition to other constraints lies in the region below 80 GeV. An 
analysis of the maximum and the minimum neutralino-proton scalar cross section for the 
allowed parameter space including the effect of a new determination of the pion-nucleon 
sigma term is also given. It is found that the full parameter space for this class of models 
can be explored in the next generation of proposed dark matter detectors. 
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1 Introduction 



Recently supersymmetric dark matter has come under a great deal of scrutiny due to 
the fact that the neutralino-proton cross sections for a wide class of supersymmetric 
models fall within the range that is accessible to the current and planned dark matter 
experiments 0, ^ ^. Thus some recent studies have included a variety of effects in 
the predictions of relic densities and of detection rates in the direct and in the indirect 
detection of dark matter 0. These include the effects of nonuniversality of the scalar 
masses at the unification scale|0, effects of CP violation with EDM constraints [Q, effects 
of coannihilation0, the effects of the — 2 constraint, as well as the effect of variations of 
the WIMP velocity []TU|, |TT|, , and the effects of rotation of the galaxy |TB| in the prediction 
of detection rates for the direct and the indirect detection of dark matter. In this work 
we focus on the effects of constraints of Yukawa unification on dark matter. This topic 
has largely not been addressed in the literature, except for the work of Ref. |T4| which, 
however, does not take account of gaugino mass nonuniversalities which is an important 
element of the present work. We focus on models where Yukawa unification occurs for ^ 



positive (we use the sign convention of Ref.|T^) consistent with the 6 s + 7 and the 
— 2 constraint from Brookhaven. The outline of the rest of the paper is as follows: In 
Sec. 2 we discuss the framework of the analysis. In Sec. 3 we discuss the (7^ — 2 constraint 
which affects dark matter analyses. This constraint requires a revision because of a recent 
reevaluation of the light by light hadronic contribution to — 2. In Sec. 4 we discuss the 
results for the satisfaction of the relic density limits under Yukawa unification constraint 
with yU > 0. In Sec. 5 we discuss the neutralino-proton cross sections including the effect 
of a new determination of the pion-nucleon sigma term. In Sec. 6 we give conclusions. 



2 Theoretical Framework 

The primary quantity of interest in the study of dark matter is ^x^"^ where — Pxl Pci 
where is the neutralino relic density, pc = SH^/SttGn is the critical matter density, 
and h is the value of the Hubble parameter Hq in units of 100 km/sMpc. Experimentally 
the limit on h from the Hubble Space Telescope is h = 0.71 ± 0.03 ± 0.07|[T^. The total 
Q = Qm + where Qm is the total matter density and is the dark energy density. 
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For we assume the simple model fim = ^b + ^x^ where fi^ is the baryonic component 
and is the neutralino component which we assume constitutes the entire dark matter. 
Using the recent analysis of Qm which gives |ll7[| Qm = 0.3 ±0.08 and assuming Qb — 0.05, 



one finds 

fi^^/i^ = 0.126 ±0.043 (1) 

With the above numerics and using a rather liberal error corridor we have the following 
limits on Q-^h'^ 

0.02 < Q^h'^ < 0.3 (2) 

In the determination of the neutralino relic density we use the standard techniques and 
compute Q^h"^ using the formula 

/T \ ^ / T \ 3 /v^/^ 



" "'^^ ^ T^J V2.73; J{xf) 

rxf 



, 3 



J (xf) = n dx {av) (x) GeV'2 (3) 

•J 

In the above (^t^J is the reheating factor, Nf is the number of degrees of freedom at the 
freeze-out temperature Tj and Xf = kTf/m^. The determination of J(a;/) is carried out 
using the accurate techniques developed in Ref . |I8| . 

It is known that gaugino mass nonuniversalities can significantly affect neutralino relic 
densities and dark matter searches. Specifically in Ref. an analysis of the effects of 
nonuniversalities of the gaugino masses on dark matter was carried out in the framework 
of SU(5) grand unification and in D brane models. We note in passing that there is no 
rigid relationship between the ratios of SU{3) x SU{2) x f/(l) gauge coupling constants 
at the GUT scale (Mg ~ 2 x lO^^ GeV) and the ratio of SU{3), SU{2), and f/(l) 
gaugino masses at the GUT scale. The ratios of the gauge coupling constants at the 
GUT scale are determined purely by the GUT group, while the gaugino masses are soft 
SUSY breaking parameters which involve both GUT and Planck scale physics. This 
topic has been discussed in several works (see for example, Refs. JTU|, pOj , |2T[ and the 
references therein) . For the present analysis we assume nonuniversality of gaugino masses 
and impose unification of gauge coupling constants at the GUT scale. Returning to the 
general structure of the gaugino masses one finds that for the case of SU(5) the gaugino 
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mass terms can arise from any of the representations that he in the symmetric product 
of 24 X 24. Since 



(24 X 24),j,„ = 1 + 24 + 75 + 200 (4) 

one finds that in general the gaugino masses are nonuniversal at the GUT scale arising 
from nonuniversalities due to the 24, 75 and 200 plet on the right hand side and one may 
in general write the SU{3) x SU{2) x U{1) gaugino masses as a sum 

rhi{MG) =m^Y^ Crul (5) 

r 

where are characteristic of the representation r and Cr are relative weights of the 
representations in the sum. Specifically the SU{3), SU{2), and f/(l) gaugino masses at 
the GUT scale for different representations have the following ratios||20|: M3 : M2 : Mi 
are in the ratio 2 : —3 : —1 for the 24 plet, in the ratio 1:3: —5 for the 75 plet, and in 
the ratio 1 : 2 : 10 for the 200 plet. The 24 plet case is of special interest for reasons which 
we now explain. It is known that the sign of the supersymmetric contribution to is 



directly correlated with the sign of in mSUGRA[23| and in other models which share 
the same generic features as mSUGRA. Thus in mSUGRA one finds that a^^^^ > for 
yU > and a^J^^^ < for fi < 0. Since experiment indicates aj^^*^^ > one infers that the 



sign of fi is positive|2^. The fi sign is also of relevance for the satisfaction of the b —>■ s + j 



constraint. It is known that the 6 — > s + 7 constraint favors a positive value of [2B 
However, a problem arises because b — t unification appears to favor a negative value of 



/ipTl, |28| . This is so because the supersymmetric correction to the b quark mass from 
the dominant gluino exchange contribution [^] is negative for /i negative and a negative 
contribution to the b quark mass helps 6 — r unification. Several works have appeared 



recently which explore ways to help resolve this problem [pOj pl| , |3^ , p3i . Specifically it 



was pointed out in Ref.|^3[ that gaugino mass nonuniversalities possess a mechanism 
which can generate a negative contribution to the b quark mass for a positive fi. In the 
context of SU(5) the mechanism arises from the gaugino mass ratios coming from the 24 
plet of SU(5) in Eq.(4). For this case it was shown that an a^^*^^ > naturally leads 
to a negative correction to the b quark mass even for /i > 0. This phenomenon comes 
about because the gluino exchange contribution to the b quark mass is proportional to 
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/im3 and the opposite sign correlation between m2 and rh^ natually leads to a negative 
contribution to the b quark mass. In this case one finds that all the constraints including 
b — T unification and the 6 — > s + 7 constraint are easily satisfied. 

One can investigate the phenomenon discussed above also in the context of SO(IO). 
Here in general the gaugino masses will lie in the symmetric SO (10) irreducible represen- 
tations of 45 X 45 where 

(45 X 45),y„ = 1 + 54 + 210 + 770 (6) 

In this case the nonuniversalities of the gaugino masses arise due to the 54, 210 and 770 
plets on the right hand side. Specifically, here one finds that the 54 plet case can give rise 
to patterns of SU{3), SU{2), and f/(l) gaugino masses which look interesting for b — t — r 
unification. Thus for the symmetry breaking pattern 5*0(10) —>■ SU{A) x SU{2) x SU{2) 
SU{3) X S{2) X f/(l) one finds that the SU{3), SU{2), U{1) gaugino masses at the 
GUT scale are in the ratio M3 : M2 : Mi = 1 : -3/2 : -ig. The SU{3) and SU{2) 
gaugino masses here have opposite signs and are similar to the 24 plet case. Thus this 
case looks favorable for generating negative corrections to the b quark mass consistent 
with a^J-^^^ > and for the unification of Yukawa couplings for /i > 0. There is another 
pattern of symmetry breaking which also looks interesting. Thus the symmetry breaking 
pattern 50(10) ^ SU{2) x 50(7) ^ 5f/(3) x 5(2) x f/(l) yields 5f/(3), 5f/(2), and 
f/(l) gaugino masses at the GUT scale in the ratio M3 : M2 : Mi = 1 : —7/3 : 1|5^. Here 
again the SU{3) and SU{2) gaugino masses have opposite signs and it appears possible 
to get negative corrections to the b quark mass necessary for b — t unification consistent 
with a^^^^ > and /i > 0. Thus we will discuss the following three cases in this paper: 
(a) 24 plet case: This is the model where the nonuniversalities originate from the 24 plet 
of SU(5) where M3 : M2 : Mi = 2 : —3 : —1; (b) 54 plet case: This is the model where 
the nonuniversalities originate from the 54 plet of SO(IO) and the symmetry breaking 
pattern gives M3 : M2 : Mi = 2 : —3 : —1; (c) 54' plet case: This is the model where the 
nonuniversalities originate from the 54 plet of SO (10) and the symmetry breaking pattern 
gives M3 : M2 : Mi = 2 : —7/3 : 1. The question that remains to be explored is what 
happens to supersymmetric dark matter in these models in the parameter space which is 
consistent with Yukawa unification and consistent with other constraints. We will discuss 
this topic in Sees. 4 and 5 after reviewing the revised constraint on — 2 in Sec. 3. 
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3 The Revised — 2 Constraint 

The recent Brookhaven experimental result gives[^] a^^*' = 11659203(15) x 10^^° where 



O'li = {9(1 — 2)/2. The standard model prediction for this quantity consists of|^ the 
0{a^) qed correction, the one and the two loop electro- weak corrections and the hadronic 
correction!^. The hadronic correction has been rather controversial ||38|] . It consists of the 



O(a^) and 0{a^) hadronic vacuum polarization, and the light-by-light hadronic contribu- 
tion. For the light by light hadronic contribution two previous analyses gave the following 
values: aj^^y^''"''^ = -7.9(1.5) x 10-i°[|9l and aj^^y^'^"'^ = -9.2(3.2) x 10-^°[|g. These give 



an average of a^f'^^'^"''^ = —8.5(2.5) x 10~^°. Using the O(a^) and 0{a^) hadronic vacuum 
polarization analysis of Ref.|5^ and the average light by light contribution as discussed 
above one finds a|f^^ — a^^^ = 43(16) x 10~^°. However, a very recent reevaluation of light 
by hght contribution finds |^ ^LbyL-Md ^ +8.3(1.2) x 10^^° which although essentially 



of the same magnitude is opposite in sign to the previous determinations. Spurred by the 
above, Hayakawa and Kinoshita^2[ reexamined the light by light contribution and found 



an error in sign in the treatment of the e— tensor in the algebraic manipulation program 
FORM used in their analysis. Their revised value of a^^'^'"-'^ = +8.9(1.54) x 10~^° is now 



in good agreement with the analysis of Ref.[41]. The average of the two evaluations gives 



^LbyL;had — _|_g gj']^^ x 10 Correcting for the above one finds 



^exp _ ^SM ^ 26(16) X 10-^° (7) 

which is a 1.6(7 deviation between experiment and theory. We discuss now the implications 
of this constraint relative to the case when one had a^^^^ — a^*^ = 43(16) x 10"^''. For the 
case when the a'^^^ — aj^*^ difference was taken to be 43(16) x 10~^° one found using a 2a 
error corridor interesting upper limits on the soft SUSY parameters and specifically for 
the mSUGRA case one found that the upper limits on mo and mi were mo < 1.5 TeV 

2 

and mi < 800 GeV for a range of tan/3 values of tan (3 < 55||2^. These ranges are well 
within the discovery limit of the Large Hadron Collider [^. 

Since the reevaluated difference af,^^ — a^,^^ is now less than 2cr, we consider a reduced 
error corridor to obtain meaningful constraints. We give here an analysis under two 
separate assumptions for the error corridor: one for 1.5cr and the other for la. The 



results of the analysis with these error corridors are exhibited in Figs. 1(a)- 1(e) for tan/3 
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values of 5, 10, 30, 45 and 55 and /i > 0. The 1.5(7 case of Fig. |l(a)| with tan/? = 5 gives 
the upper hmits tuq < 850 GeV and mi < 800 GeV. However, here the lower limit of 
Higgs boson mass indicated by the LEP data lies outside of the allowed parameter space. 
The 1.5(7 case of Fig. 1(b) with tan/5 = 10 gives the upper limits tuq < 1300 GeV and 
mi < 1100 GeV. Here the parameter space includes the lower limit of the Higgs boson 

2 

mass indicated by the LEP data. For the case of Fig. |l(c)| with tan j3 = 30 one finds the 
upper limit of mo < 2500 GeV for the 1.5a case which is on the borderline of the reach 



of the LHC||4^ and most likely beyond its reach. However, for the la case one finds the 
upper limits mo < 1000 GeV and mi < 800 GeV which lie well within the discovery 
potential of the LHC. A similar situation holds for tan/5 = 45 and for tan/? = 55. For 
the tan/5 = 45 case of Fig. p.(d)| the upper limit for mo is rriQ < 2700 GeV for the 1.5a 
case while the la case gives mo < 1300 GeV and mi < 825 GeV, which lie well within 
the reach of the LHC. Similarly the tan /? = 55 plot of Fig. |l(e)| gives mo < 2500 GeV for 
the 1.5(7 case but the la case gives mo < 1450 GeV and mi < 625 GeV. Again while the 

2 

1.5cr upper limits are on the border line of the reach of the LHC and most likely beyond 
its reach, the la upper limits lie well within reach of the LHC. Since the upper limits for 
the 1.5(7 case appear to cross the usual naturalness limits (see, e.g., Ref.|Q) at least for 
values of tan/? > 30 we impose the constraint of la error corridor around the mean for 
a^^P — a^^^ as given by Eq.(7) for the analysis of Sees. 4 an 5. The upper limits implied by 
the latter constraint lie well within the naturalness limits. 



4 Yukawa Unification and Relic Density Analysis 

We turn now to the main theme of the paper which is the analysis of dark matter under the 
Yukawa unification constraints. As discussed in Ref. p3| we define the Yukawa coupling 
unification parameter Sij for the Yukawa couplings Aj an Xj so that 

where Aj-,- = (Aj + Xj)/2 and 6ij defines the degree of unification. As is well known dark 
matter analysis is very sensitive to the 6 s + 7 constraint. There are several recent 
experimental determinations of 6 — s + 7 i.e., CLEO[^], BELLE[^ and and we 
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use their weighted mean. Analyses of the theoretical prediction of the Standard Model 
including the leading order and the next to leading order corrections for this branching 
ratio have been given by several authors In our analysis we use a 2a corridor in the 
difference between experiment and the prediction of the Standard Model to constrain our 
theoretical analysis of the supersymmetric contribution. First we discuss the SU(5) case 
where we consider the gaugino mass nonuniversality at the GUT scale from the 24 plet 
representation as discussed in the paragraph following Eq.(5). The remaining soft SUSY 
breaking parameters in the theory are assumed universal. 

We begin by discussing the allowed parameter space in the mo — C24 *mi plane under 
the Qfj. — 2 constraint. The results are exhibited in Fig. P(a)| for values of tan/? of 5, 10, 
30 and 40. The top gray regions correspond to disallowed areas via radiative electroweak 
symmetry breaking requirement. The bottom patterned regions in Fig. 2(a) for the cases 
of tan/5 = 5,10 and 30 are typically eliminated via the stability requirement of the Higgs 
potential at the GUT scale. Part of the region with large |c24mi/2| and large mo bordering 
the allowed (white) region for tan/3 = 30 is eliminated via the limitation of the CP-odd 
Higgs boson mass turning tachyonic at the tree level. For tan/3 = 40 most of the region 
(patterned and shaded) to the right of the allowed white small region is eliminated because 
of Afc going to the non-perturbative domain due to a large supersymmetric correction 
to the bottom mass. Regarding the SO(IO) well known one needs to use 

nonuniversality of the Higgs boson masses at the GUT scale to achieve radiative breaking 
of the electroweak symmetry. For the analysis here we use the nonuniversal Higgs scalar 
masses so that m|^^ = l.Smg and — O.Smg. The result of analysis of the allowed 
parameter space in the 54 plet case is given in Fig. P(b)| for values of tan/3 of 5, 10, 30 
and 45. The regions with patterns are discarded for reasons similar to as in Fig. ^(a) 



There is no discernible change in these results due to modest variations (up to 50%) in 
the assumed nonuniversality (i.e., deviations of m^^ and m^^^ at Mq from mg) of the 
soft parameters in the Higgs boson sector needed to accomplish radiative breaking of the 
electroweak symmetry. 

We give now the relic density analysis. As a guide we use the unification criterion 



6bT < 0.3. In Fig. |3(a)| we plot f^^^^ tan/3 for the following range of parameters: 

< mo < 2000 GeV, -1000 GeV < C24mi < 1000 GeV, -6000 GeV < Ao < 6000 GeV 

2 

and fi > 0. The small crosses are the points that satisfy the — 2 constraint, the filled 
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squares additionally satisfy the 6 — > s + 7 constraint and the filled ovals satisfy all the 
constraints, ie., the — 2 constraint, the 6 — s + 7 constraint and the h — t unification 
constraint with 5bT < 0.3. One finds that there exist significant regions of the parameter 
space as given by filled ovals where all the constraints are satisfied. The horizontal lines 
indicate the allowed corridor for the relic density as given by Eq.(2). A plot of f^^^^ 



tuq for exactly the same ranges of the parameter space as in Fig. |3(a)| is given in Fig. ^(b) 



A similar plot of ^^^'^ ^ function of C24 * mi is given in Fig. g^. In Fig. gd] 



we 



give a plot of ^^^^^ ^0 and in Fig. |3(e)| we give a plot of Q^h"^ as a function of the 
LSP mass m^. The paucity of points in the region around the neutralino mass of 45 GeV 
in the allowed corridor of relic density in Fig. |3(e)| is due to the rapid s-channel Z pole 
annihilation and also due to the s-channel Higgs pole annihilation in the region below 
rriy, ~ 60 GeV. Finally in Fig. 3(f) we give a plot of ^^^"^ ^br- One finds that there 
exist regions of the parameter space where b — t unification even at the level of a few 
percent consistent with the relic density and other constraints can be satisfied. 

We discuss next the SO(IO) case with gaugino mass nonuniversality of the type M3 : 
M2 : Ml = 1 : — I : — 1 as given in the paragraph following Eq.(6). As noted earlier 
in this case the pattern of relative signs of the gaugino masses is similar to that for the 
24 plet case. In the analysis we impose not only the b — t Yukawa unification but also 
b — t and t — t Yukawa unification Extrapolating from the SU(5) case we impose the 
following constraints on 6i,r, 5u and dtr'- all 5ij < 0.3. In Figs. |4(a)| we plot Q^h"^ vs tan/3 
for for the same range of parameters as in Fig. |3(a)| . The symbols used in Fig. [4(a) 



I.e., 



the small crosses, the filled squares and the filled ovals also have the same meaning as in 
Fig. |3(a)| except that the filled ovals now include all the Yukawa unification constraints, 
i.e., 6bT, ^bt, ^tr < 0.3. Fig. [4(a)| shows that there exist significant regions of the parameter 
space as given by filled ovals where all the constraints including the Yukawa unification 
constraints are satisfied. The horizontal lines indicate the allowed corridor for the relic 
density as given by Eq.(2). In Fig. 4(b) we give a plot of ^^^"^ "^0 for exactly the 
same ranges of the parameter space as in Fig. [4(a)| . A similar plot of ^^^"^ ^ function 
of C54 * mi is given in Fig. [4(c)| . In Fig. [4(d)| we give a plot of fi^^/i^ vs Aq and in 



Fig. [4(e)| we give a plot of ^ function of the LSP mass m^^. Again the paucity 

of points in the region around the neutralino mass of 45 GeV in the allowed corridor of 
relic density in Fig. [4(e)| is due to the rapid Z pole annihilation. Finally in Fig. [4(f)| we 
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give a plot of ^x^"^ ^br and similar plots exist for Q^^'^ vs 6bt and vs 6tT but 

are not exhibited. Interestingly in this case one finds a high density of points where the 
relic density constraint consistent with Yukawa unification at the level of a few percent is 
satisfied. An analysis of the gaugino-higgsino content of the neutralino over the parameter 
space of the model consistent with b — t unification and all the other constraints can be 
gotten by examining the expansion of the LSP so that 

X = aB + pW3 + iHi + 6H2 (9) 

For the 24 plet case one finds that typically over most of the parameter space > 0.75 

while for the 54 plet case over most of the parameter space one has + (3"^ > 0.03. While 
the dominant component in all cases is the Bino|Q, in the 24 plet case one could also 



have a significant higgsino component to the LSP while for the 54 plet case the Bino 
purity of the LSP is rather high. We have carried a similar analysis for the relic density 
for the second SO (10) case discussed in Sec. 2 where M3 : M2 : Mi = 1 : — | : 1 (we label it 
as the 54' case). We do not exhibit the details as this case appears marginal in the sense 
that the allowed neutralino mass range is very small after the LEP limit on m^^ > 32.3 
GeV is imposed. One should keep in mind, however, that the LEP limit is a generic limit 
and is not deduced specifically for the model under discussion. Still this case has rather 
low neutralino mass upper limit in any case. For the sake of completeness we will discuss 
the neutralino-proton cross sections for this case also in Sec. 5. 

The allowed mass ranges for the three cases discussed above, i.e., 24, 54 and 54' cases 
are given in Table 1. The spectrum of Table 1 satisfies all the desired constraints, i.e., 
— 2, 6 — > s + 7 and Yukawa unification constraints as discussed above. A remarkable 
aspect of Table 1 is that the b — r unification constraint implies a rather low Higgs boson 
mass. Since tan jS values for scenarios with b — t unification imply rather high values of 
tan/5, the experimental lower limits from LEP for large tan/5 are rather low, i.e., rrih > 91 
GeVfS^. Thus the higgs mass ranges listed in Table 1 are all consistent with the current 



experimental limits. Further, these limits are tantalizingly close to observation at RUNII 
of the Tevatron. Further, the spectrum for all the three cases listed Table I is accessible to 
the LHC As noted earlier the 54' case has only a very narrow allowed range in neutralino 
mass and could be tested or eliminated by data with a modest improvement in energy. 
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Table 1: Sparticlc mass ranges for 24, 54, and 54' cases 



Particle 


24 (GeV) 


54 (GeV) 


54' (GeV) 




32.3 - 75.2 


32.3 - 81.0 


32.3 - 33.4 


xl 


96.7 - 422.5 


94.7 - 240.8 


145.7 - 153.9 


xl 


110.5 - 564.3 


301.5 - 757.1 


420.9 - 633.8 


xl 


259.2 - 575.9 


311.5 - 759.7 


427.6 - 636.9 


Xl 


86.9 - 122.6 


91.6 - 210.8 


1 15.8 - 153.9 


xt 


259.9 - 577.2 


315.1 - 761.6 


430.7 - 639.2 


9 


479.5 - 1077.2 


232.5 - 580.3 


229.8 - 237.4 




299.7 - 1295.9 


480.5 - 1536.8 


813.1 - 1196.3 




355.1 - 1309.3 


489.8 - 1482.7 


835.3 - 1237.6 




203.5 - 1045.1 


294.2 - 1172.6 


579.4 - 863.7 




349.6 - 1180.9 


422.6 - 1311.7 


704.6 - 1018.3 


Ul 


533.6 - 1407.2 


566.7 - 1506.4 


822.9 - 1199.8 


U2 


561.1 - 1443.0 


584.7 - 1544.6 


849.6 - 1232.6 


di 


535.1 - 1407.5 


580.3 - 1546.2 


845.1 - 1232.5 


d2 


566.7 - 1445.2 


590.1 - 1546.7 


853.3 - 1235.2 


ii 


369.9 - 975.2 


271.5 - 999.6 


513.7-819.9 


h 


513.7- 1167.6 


429.4 - 1107.4 


599.4 - 848.2 


hi 


488.2 - 1152.8 


158.1 - 1042.0 


453.2 - 749.9 


h 


532.3 - 1207.0 


396.6 - 1159.2 


610.5 - 880.4 


h 


104.3 - 114.3 


103.8 - 113.3 


108.1 - 110.9 


H 


111.9 - 798.8 


151.5 - 1227.6 


473.4 - 831.9 


A 


110.5 - 798.8 


151.4 - 1227.6 


473.4 - 831.9 


1^ 


96.0 - 559.5 


291.1 - 752.7 


413.1 - 628.4 



5 Maximum and minimum neutralino-proton cross 
section 

In the analysis of neutralino-proton cross section a^^p we restrict ourselves to the CP 
conserving case. Here the x~P scattering is governed by the four Fermi interaction Ceff = 
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Xli^l^Xqi^iAPi + BPR)q + Cxxmgqq + i3>X75X"^g<?75g + Exh^X^ggq + Fxxmqqh^q- We 
are specifically interested in neutralino scattering from heavy targets. This scattering is 
dominated by the scalar interactions and in this case the x ~ P cross-section is given by 

V, ^ „„„ 2 

7T 



a^,{scalar) = ^{ J] /fC. + ^(1 - ^ /f) ^ C^)' (10) 



i=u,d,s i=u,d,c a=c,b,t 

In the above /f are the (i=u,d,s) quark densities which are defined by m^/f =< p\mgiqiqi\p >, 
and fir is the reduced mass. The scalar interaction parametrized by C arises from several 
sources: from s channel contributions from the higgs h^, exchanges and from t channel 
contributions from the sfermion exchange so that C = Cho + Cj^o + Cj. It was shown in 
Ref. |I3 that it is convenient to parameterize the form factors //^'"^ such that 



f{u,d) = ^{u,d) i^n + md) ^(1 ± O^^TrNmp^, /f 



ms{mu + md) ^(1 - a;)cr^Armp ^ 



{u,d) 



m 



(u,d) 



{rriu + md) ^{1 T O^ttnTtiJ-, = ms{mu + nid) ^{l - x)a^Nm^ (11) 



where a^^jq^ x and ^ are defined by cr^^=< p\2^^ {mu+rnd){uu+dd\p >, ^ =< p\uu — dd\p > 
(< p\uu + dd\p >)^^, and x = (yQ/(y-nN=< p\uu + dd — 2ss\p > {< p\uu + dd\p >)^^ . Quark 
densities for the neutron are related to the proton quark densities by|jl9|/^/| = fufd- 



Baryon mass splittings can be used to determine the ratio ^/x and one finds|T^ ^/x = 
0.196. Using various determinations of (Tq and (Jt^n-, x was estimated in Ref.|jl9|] to be 
X = 0.67 ± 0.18 which gives ||19|| ^ = 0.132 ± 0.035. Using the current data on the quark 
masses one finds fP = 0.021 ± 0.004, = 0.029 ± 0.006, and /f = 0.21 ± 0.12 and 
= 0.016 ± 0.003, f2 = 0.037 ± 0.007, and = 0.21 ± 0.12. 

It has been pointed out recently]^ that an analysis of a7rAr[§3| using new pion-nucleon 
scattering data[Q leads to a significantly larger neutralino-nucleon cross section. Thus 
the new determination of o"7riv[Q lies in the range 55 MeV < a-n-N < 73 MeV which 
is much larger than the previous determinations (see, e.g., Ref. ||19||). Using the new 



determination of a-j^M and repeating the analysis of Ref.|T^] we find x = 0.55 ± 0.12 and 
C, = 0.108 ±0.024. These lead to the following new determinations for the quark densities 



/„^ = 0.027 ±0.005, fd 
/: = 0.022 ±0.004, 



0.038 ±0.006, /f 
0.049 ±0.007, /; 



0.37±0.11 
0.37 ±0.11 



(12) 



We use these new quark densities in our numerical analysis. In Fig. |5(a)| we exhibit the 
neutralino-proton scalar cross section vs the neutralino mass for the nonuniversal case 
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of Fig. 1^. The gaps in Fig. |5(a)| are due to the rehc density constraint as can be seen 
from Fig. |3^ and Fig. [3(e]i The DAMA regionfl]], the lower hmit achieved by CDMSg] 
and the future lower limits that may be achieved @, |^ are also exhibited. First, one finds 
that the allowed neutralino range is significantly reduced in this scenario with mass range 
limited to less than 65 GeV. Second, one finds that the parameter space of the model 
can be fully probed by the proposed future dark matter detectors[^, ^. In this model 
a^p{scalar) lies in the range 



4 X 10"^^ 



cm) 



< a^p{scalar) < 4 x 10 



-41 



cm) 



(13) 



In Fig. |5(b)| we give a plot of a^p{scalar) vs tan/3 which shows that the upper limits of 
a^p{scalar) are strongly dependent on tan/3 as expected. An analysis of a^p{scalar) vs 
m^ for the nonuniversal SO(IO) gaugino mass case of Fig. |4(c)| is given in Fig. |5(c)| . Here 



the neutralino mass range extends up to 80 GeV. As for the case of Fig. 5(a) the gaps in 
Fig. |5(c)| are due to the relic density constraint as can be seen from Fig. P(c)| and Fig. |3(e)| . 
Again as in the 24 plet case all of the parameter space of this model can be fully probed 
by the proposed future dark matter detectors|Q, ^. In this model a^p{scalar) lies in the 
range 



7 X 10"^^ 



cm) 



< a^p{scalar) < 1 x 10 



-41 



cm) 



(14) 



In Fig. |5(d)| we give a plot of a^p{scalar) vs tan (3. Here since tan (3 does not vary over a 
wide range one does not see a large enhancement of a^p{scalar) with tan/? in this limited 
range. Finally, we discuss a^p{scalar) for the 54' case. In Fig. |5(e)| we exhibit (y^p{scalar) 
as a function of m^. In Fig. [5(f)| we give a plot of a^p{scalar) as a function of tan/?. The 
lower and upper limits on the scalar cross section in this case are very similar to the 54 
plet case of Figs. |5(c)| and ^(d)| . However, as discussed in Sec. 4 imposition of the lower 
limit of 32.3 GeV on the neutralino mass eliminates most of the parameter space of this 
model. 



6 Conclusions 

In this paper we have given an analysis of supersymmetric dark matter under the con- 
straint of Yukawa coupling unification with /i > 0. The constraints of 6 — s + 7 and the 
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revised — 2 constraint taking account of the recent reevaluation of the hght by hght 
hadronic correction were also imposed. The analysis was done exploiting the recently 
discovered region of the parameter space which utilizes nonuniversal gaugino masses and 
leads to negative corrections to the b quark mass necessary for Yukawa coupling unifica- 
tion with /i > 0. We considered scenarios with SU(5) and SO(IO) unifications. Within 
SU(5) we considered nonuniversalities arising from the 24 plet representation of SU(5) 
which allow for significant regions of the parameter space consistent with h — r unification 
so that 5},r < 0.3 for /i > consistent with other constraints. This scenario limits the 
neutralino mass range to lie below 65 GeV and within this range a significant part of 
the parameter space is consistent with the relic density constraint. An analysis of the 
neutralino-proton scalar cross section reveals that the allowed range of cross sections can 
be fully probed by the proposed future dark matter detectors. Within SO (10) we consid- 
ered nonuniversalities arising from the 54 plet representation of SO(IO) which allow for 
significant regions of the parameter space consistent with for, ht and tr unification con- 
straints such that 5i,r,5}yt,5tr < 0.3 for /i > consistent with other constraints. In this 
case one finds that the neutralino mass range extends till 80 GeV and again the analysis 
of neutralino-proton scalar cross section shows that the allowed range of cross sections 
can be fully probed by the proposed future dark matter detectors. One of the important 
features of models with h — r unification explored here is a relatively low lying light Higgs 
boson with mass lying below 115 GeV. This mass range would certainly be explored by 
RUNII of the Tevatron. Further, the entire sparticle spectrum predicted in the class of 
models with h — r unification discussed here would be accessible at the LHC. It would be 



interesting to explore SUSY signals such as the trileptonic signal||55| at colliders from this 
model. But such an investigation is beyond the scope of this paper. 
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Figure Captions 
Fig. I: 

Fig. |l(a)| : Allowed aj^^*^^ regions corresponding to the 1.5a and the Icr constraints for 



tan/3 = 5. Similar analyses for tan/? = 10,30,45,55 are given in Fig. 1(b) , Fig. l(c 



Fig. p.(d)| and Fig. p.(e)| respectively. The top left gray regions do not satisfy the radia- 
tive electroweak symmetry breaking requirement whereas the bottom patterned regions 
are typically discarded by stau becoming the LSP. For large tan (3, 45 or 55, the bottom 
patterned region near the higher mi/2 side and on the border of the white allowed regions 
are discarded because of CP-odd Higgs boson turning tachyonic at the tree level. 



Fig. 0: 

Fig. |2(a)| : Allowed — 2 regions corresponding to 1.5a and la constraints for nonuniver- 
sal gaugino mass scenario of the SU(5) 24 plet case. A discussion of the discarded regions 
in the top and the bottom parts of the figures is given in the text in Sec. 4. 
Fig. ^(b)| : Allowed — 2 regions corresponding to 1.5cr and la constraints for nonuni- 
versal gaugino mass scenario of SU(IO) 54 plet case. Here the nonuniversal Higgs scalar 
parameters are given by m|^^ = 1.5mQ and m^^^ = 0.5mQ. A discussion of the discarded 
regions in the top and the bottom parts of the figures in given in Sec. 4. 



Fig: |: 

Fig. |3(a)| : Plot of ^x^'^ tan/5 for the SU(5) 24-plet case with the inputs < rriQ < 2 
TeV, —1 TeV < C24mi/2 < 1 TeV, —6 TeV < Aq < 6 TeV and /i > 0. The small crosses 
satisfy the — 2 constraints, the (blue) filled squares additionally satisfy the 6 ^ s + 7 
limits and the (red) filled ovals satisfy all the constraints, i.e., the g^ — 2 constraint, the 
6 ^ s + 7 constraint, and b — r unification at the level Sbr < 0.3. The two horizontal lines 
refer to the limits of Eq. (2) 

Fig. 3(b) : Plot of Q^h"^ vs mo with all the same parameters as in Fig. 3(a) and with 



tan P < 55. Symbols have the same meaning as in Fig. |3(a) 



Fig. |3(c)| : Plot of vs C24 * mi with all other parameters the same as in Figs. |3(a) 



and |3(b)| . Symbols have same meaning as in Fig. |3(a' 



Fig. |3(d)| : Plot of ^^^"^ ^0 with all other parameters the same as in Figs. P(a)| , |3(b) 
and |3(c)| . Symbols have the same meaning as in Fig. |3(a 
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Fig. P(e)|: Plot of Q^h'^ vs the LSP mass with all other parameters the same as in 



Figs. p(a)| to |3(d)| . The small plus symbols refer to valid parameter points with no con- 
straints, and (red) filled ovals refer to satisfying all the constraints i.e., the constraints 
from (yf^ — 2, 6 — i> s + 7 and 6br < 0.3. 

Fig. P(f )| : Plot of ^x^"^ ^br with all other parameters the same as in Figs. |3(a)| to P(d) 
Symbols have the same meaning as in Fig. |3(a) 



Fig: |: 



Fig. [4(a)| : Plot of r^^^/i^ vs tan/5 for the SO(IO) 54-plet case with inputs < < 2 
TeV, —1 TeV < 0^41711/2 < 1 TeV, —6 TeV < Aq < 6 TeV and /i > 0. The nonuniversal 
Higgs scalar parameters are given by m|^^ = LSttIq and m^^^ = 0.5mQ. The small crosses 
satisfy the — 2 constraints, the (blue) filled squares additionally satisfy the 6 — s + 7 
limits and the (red) filled ovals satisfy all the constraints, i.e., the — 2 constraint, the 
6 s + 7 constraint, and unification of Yukawa couplings so that 6bT, Su, ^tr < 0.3. The 
two horizontal lines refer to the limits of Eq.(2) 

Fig. |4(b)| : Plot of VL^h"^ vs itlq with all the same parameters as in Fig. [4(a)| and with 



tan (3 < 55. Symbols have the same meaning as in Fig. [4(a 



Fig. |4(c)| : Plot of Q^^"^ ^^54 * mi with all other parameters the same as in Figs. [4(a^ 
and 4(b) . Symbols have the same meaning as in Fig. 4(a^ 



Fig. 4(d)| : Plot of ^x^"^ ^0 with all other parameters the same as in Figs. Ei], gb] 
and [4(c)| . Symbols have the same meaning as in Fig. |4(a 



Fig. [4(e)| : Plot of Q^h vs the LSP mass with all other parameters the same as in 
Figs. [4(a)| to [4(d)| . The small plus symbols refer valid parameter points with no con- 
straints, and (red) filled ovals refer to satisfying all the constraints i.e., the constraints 
from (yf^ — 2, 6 — i> s + 7 and 6bT, SbAr < 0.3. 

Fig. [4(f )| : Plot of ^x^"^ ^br with all other parameters the same as in Figs, ga]] to gd] 
Symbols have the same meaning as in Fig. 4(a) 



Fig. §: 



Fig. ^(a)| : Plot of the neutralino-proton scalar cross section cr^p vs the lightest neutralino 
mass rrix for the SU(5) 24 plet case with the range of the parameters given in Figs. |3(a 
to P(d)| satisfying all the desired constraints including the b — t unification constraint so 
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that 5bT < 0.3. The small crosses satisfy the — 2 constraints, the (blue) filled squares 
additionally satisfy the b s+7 limits and the (red) filled ovals satisfy all the constraints, 
i.e., the — 2 constraint, the 6 — > s + 7 constraint, and 6bT < 0.3. The area enclosed by 
solid lines is excluded by the DAMA experiment [Q , the dashed line is the lower limit from 
the CDMS experiment]^], the dot-dashed line is the lower limit achievable by CDMS in 
the future and the dotted line is the lower limit expected from the proposed GENIUS 
experiment Q . 

Fig. |5(b)| : Plot of the neutralino-proton scalar cross section cr^p vs tan/3 for the SU(5) 24 
plet case with the same range of parameters as given in Figs. |3(a)| - |3(d)| satisfying all the 
desired constraints including the b — t unification constraint so that 6bT < 0.3. Symbols 
have the same meaning as in Fig. |5(a). 



Fig. ^(c)| : Plot of the neutralino-proton scalar cross section a^p vs the neutralino mass 
for the SO(IO) 54 plet case with the same range of parameters as given in Figs. |4(a 



4(d)| satisfying all the desired constraints including the b — t, b — t and t — t unification 
constraint so that 6bT, Sbt, ^tr < 0.3. The small crosses satisfy the g^ — 2 constraints, the 
(blue) filled squares additionally satisfy the 6—^5 + 7 limits and the (red) filled ovals 
satisfy all the constraints, i.e., the g^ — 2 constraint, the b ^ s + j constraint, and Yukawa 
unifications with 6bT, ht^tr < 0.3. 

Fig. |5(d)| : Plot of the neutralino-proton scalar cross section a^p vs tan/? for the SO(IO) 
54 plet case with the same range of parameters as given in Figs. 4(a) - [4(d) satisfying all 
the desired constraints including the & — r, 6 — t and t — t unification constraint so that 



5bT-, Sbt, StT < 0.3. Symbols have the same meaning as in Fig. 5(c^ 



Fig. p(e)| : Plot of the neutralino-proton scalar cross section a^p vs the neutralino mass 
for the SO(IO) 54' plet case with the same range of parameters as given in Figs. [4(a) 



- 4(d)[ satisfying all the desired constraints including the b — r, b — t and t — t unification 
constraint so that 6bT, Sbt, ^tr < 0.3. Symbols have the same meaning as in Fig. [5(c)[ . 
Fig. p(f )[ : Plot of the neutralino-proton scalar cross section a^p vs tan/? for the SO(IO) 
54' plet case with the same range of parameters as given in Figs. [4(a)[ - [4(d)[ satisfying all 
the desired constraints including the b — r, b — t and t — t unification constraint so that 



6bT, Sbt, StT < 0.3. Symbols have the same meaning as in Fig. 5(c' 
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